Gamma ray emission and stochastic particle acceleration in 

galaxy clusters 

G. Bmnetti*, P. Blasi^ R. Cassano* and S. Gabici** 

*INAF Istituto di Radioastronomia, via P.Gobetti 101, 40129, Bologna, Italy 
^INAF Osservatorio Astrofisico di Arcetri, Firenze, Italy and Theoretical Astrophysics, Fermi National 

Accelerator Laboratory, Batavia, USA 
** Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin, Ireland 

Abstract. FERMI (formely GLAST) will shortly provide crucial information on relativistic particles in galaxy clusters. We 
discuss non-thermal emission in the context of general calculations in which relativistic particles (protons and secondary 
electrons due to proton-proton collisions) interact with MHD turbulence generated in the cluster volume during cluster 
mergers. Diffuse cluster-scale radio emission (Radio Halos) and hard X-rays are produced during massive mergers while 
gamma ray emission, at some level, is expected to be common in galaxy clusters. 

Keywords: Particle acceleration, Turbulence, Galaxy clusters. Gamma rays. 
PACS: 96.50.Pw, 96.50.Tf, 98.65.Cw, 95.85.Pw. 



INTRODUCTION 

Clusters of galaxies contain « lO'^ Mq of hot (10^ K) 
gas, galaxies, dark matter and non-thermal components. 

The origin of non-thermal components is likely con- 
nected with the cluster formation process : a fraction of 
the energy dissipated during cluster mergers is expected 
to be channelled into the acceleration of particles via 
shocks and turbulence that lead to a complex popula- 
tion ofprimary electrons and protons in the IGM (e.g., 
ISIISIIISIISIISII [H). Theoretically relativis- 
tic protons are expected to be the dominant non-thermal 
particle components since they have long life-times and 
remain confined within galaxy clusters for an Hubble 
time (e.g. ifllll and ref. therein). Confinement enhances 
the probability to have proton-proton (p-p) collisions that 
in turns give gamma ray emission via decay of the neu- 
tral pions produced during these collisions |[T2i[T3ll . p-p 
collisions also inject secondary electrons that give syn- 
chrotron and inverse Compton (IC) emission whose rel- 
evance depends on the proton content in the IGM lfl4ll . 
The inter-galactic medium (IGM) is expected to be turbu- 
lent at some level and MHD turbulence can re-accelerate 
both primary and secondary particles via second order 
Fermi mechanisms. Turbulence is naturally generated in 
cluster mergers (e.g., |[il[il[i3,[il) and the resulting 
particle re-acceleration process should enhance the syn- 
chrotron radio emission and the IC hard X-ray emission 
by orders of magnitude. 

The energy content in the form of relativistic pro- 
tons in the IGM is still poorly constrained since present 
gamma ray observations can provide only upper Um- 
its to the gamma ray emission from galaxy clusters 



( lfl9ll : also Reimer and Perkins, this conference). On the 
other hand radio observations of Mpc-sized diffuse syn- 
chrotron emission from galaxy clusters provide crucial 
information on the relativistic electron component (e.g. 

The FERMI gamma-ray space telescope will shortly 
allow to measure (constrain) the energy content of rel- 
ativistic protons in the IGM. For this reason, starting 
from present understanding of non thermal components 
in galaxy clusters, we discuss expectations for gamma 
ray emission. 

CLUSTER-SCALE RADIO EMISSION 

In this Section we outline our starting point to model non 
thermal (including gamma ray) cluster emission. 

The most prominent examples of diffuse cluster emis- 
sion are giant Radio Halos: Mpc-scale diffuse syn- 
chrotron sources at the centre of a fraction of massive 
and merging galaxy clusters ll20l[2lll . This Mpc-scale ra- 
diation may originate from secondary electrons injected 
by collisions between relativistic and thermal protons 
in the IGM (e.g. |j22t alternatively extended ra- 

dio emission may originate from relativistic electrons 
re-accelerated in situ by various mechanisms associ- 
ated with the turbulence in massive merger events (e.g. 
nil,!!!). These two processes likely happen at the same 
time and a unified scenario that models both the injection 
and re-acceleration of secondary electrons and primary 
particles due to MHD turbulence has been investigated 
by iH. 

Radio Halos are not common: although a fairly large 



number of clusters has adequate radio follow up, they 
are presently detected only in a fraction of massive and 
merging clusters llzsi l26l EtI l28ll . Studies based on the 
analysis of present X-ray selected cluster samples with 
radio follow up allow to conclude that the fraction of 
clusters with Radio Halos depends on cluster X-ray lu- 
minositjQ, and that clusters have a bimodal behaviour in 
the 1.4 GHz radio luminosity {Pi 4) - soft X-ray lumi- 
nosity {Lx) plane with Radio Halo clusters being one or- 
der of magnitude more radio luminous than upper limits 
for clusters with no Radio Halos (|29]; Fig. 1). 

These facts suggest that Radio Halos are transient phe- 
nomena connected with cluster mergers and that some 
threshold in the mechanism for the generation of these 
sources should come into play. Unless we admit the pos- 
sibility of strong dissipation of the magnetic field in clus- 
ters, these properties cannot be easily understood in the 
case that the continuous injection of secondary electrons 
in the IGM plays the major role in the origin of these 
sources. In this case - indeed - Radio Halos should be 
long living phenomena and common, and some general 
A 4-L_t trend would be predicted for all clusters (e.g. 
iH 113, In particular, the bimodality in Fig. 1 and 
the substantial lack of clusters between the Radio Halo 
and radio quiet regions implies that the bulk of the mag- 
netic energy in clusters should be dissipated in a time- 
scale of only w 0.1 —0.2 Gyr which is challenging to 
reconcile with present understanding of cluster magnetic 
fields (e.g. |3lB. 

On the other hand, the emerging observational picture 
supports the idea that turbulent re-acceleration of rela- 
tivistic electrons may play a role in the formation of Ra- 
dio Halos in connection with cluster mergers. In partic- 
ular it is interesting to note that in the context of this 
scenario radio-quiet clusters are expected to evolve into 
Radio Halo clusters (and vice versa) in a time-scale of 
the order of the electron acceleration time-scale, «0.1- 
0.2 Gyr Spectral studies are also important : the dis- 
covery of Radio Halos with steep spectrum fsS] and the 
synchrotron cut off that is found in the spectrum of the 
Coma Radio Halo Issl [34ll imply a maximum energy in 
the spectrum of the emitting electrons at energies «GeV 
suggesting that the mechanism responsible for the accel- 
eration of electrons in the IGM is poorly efficient, con- 



' the fraction is!^ 30% for Lx > 8 ■ 10'''' erg s"' clusters and < 10% for 
clusters with Lx « 3 ■ 10'''' — 8 ■ 10'''' erg s"'; Montecarlo approaches 
show that the two fractions differ at « 4f7 level IitIi 
- this means that assuming that all massive clusters may host Radio 
Halos during their life, the Hfe-time of Radio Halos must be much 
shorther, < 1 Gyr, than the cluster life-time. 

^ even in the (worst) case of simply decaying MHD turbulence, the 
energy density of the rms field decreases with time only (about) linearly 
and > 5 eddy turnover times (Ri Gyr) are required to allow dissipation 
of the bulk of the magnetic field (Fig. 2 in t3li1 ) 
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FIGURE 1. Distribution of galaxy clusters in the Pi 4-Lx 
plane. Filled dots are galaxy clusters at z>0.2 (from GMRT- 
sample and from the literature), empty dots are clusters at lower 
redshift reported to highlight the /'1.4-Lx correlation (solid 
line). Upper limits (represented by the arrows) are GMRT clus- 
ters with no hint of cluster-scale emission and their distribution 
should be compared with that of clusters at similar redshift 
(filled dots). 



sistent with turbulent acceleration. 



EMISSION FROM GALAXY CLUSTERS 

In this Section we calculate non-thermal (multifre- 
quency) emission from galaxy clusters under the assump- 
tions that MHD turbulence plays a role in the parti- 
cle acceleration process during cluster mergers. We do 
not include the contribution to the non thermal emission 
from fast electrons accelerated at shock waves that de- 
velop during cluster mergers and accretion of matter (see 
Pfrommer, this conference). 

As already mentioned, it is believed that clusters are 
reservoir of relativistic protons that accumulate in the 
IGM during cluster life-time (e.g. |11]). Thus we model 
the re-acceleration of relativistic particles by MHD tur- 
bulence in the most simple situation in which only rela- 
tivistic protons are initially present in a turbulent IGM. 
These protons generate secondary electrons via p-p col- 
lisions and in turns secondaries (as well as protons) are 
re-accelerated by MHD turbulence. Following [24] we 
restrict to the case of Alfven wavefl and calculate the 
spectrum of particles and MHD waves and their evolu- 
tion with time by solving a set of coupled equations that 



an additional possibility is given by magnetosonic waves f33, l^ 




FIGURE 2. Broad band spectrum produced within R < I Mpc from a Coma-like cluster. Upper panels: Synchrotron (left, SZ 
decrement at high frequencies is not taken into account), and IC and n-o emission (right) calculated at t=0.5 Gyr from the injection 
of MHD turbulence in the IGM (the energy injected in Alfven modes between t=0-0.5 Gyr is ~3% of the thermal energy). Lower 
panels: Synchrotron (left) and IC and n-o emission (right) calculated at t=l Gyr after dissipation of turbulence in the IGM. In all 
panels calculations are shown assuming a ratio between the energy density of relativistic and thermal protons = 1% (dashed lines), 
0.5% (dotted lines) and 0.3% (solid lines) at t=0 (with proton spectrum 5 = 2.2) and a central cluster-magnetic field Bo = 2/iG. For 
the sake of completeness we show radio data, BeppoSAX data and EGRET upper limit for the Coma cluster (|24] and ref. therein) 
the recent VERITAS upper limit (Perkins, this meeting) and the approximate sensitivity after 1 yr of FERMI (dashed). 



give the spectrum of electrons, N^, , positrons, N^, pro- 
tons, Np, and waves. Wit '■ 
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where \dp/dt\ marks radiative (r) and Coulomb (i) 
losses, Dpp is the particle diffusion coefficient in the 
momentum space (and depends on the wave spectrum 
Wk), is the injection term of secondary leptons due to 
p-p collisions (and depends on Np), Dj^i^ is the diffusion 
coefficient in the wavenumber space, 4 is the injection 
rate-spectrum of Alfven waves at resonant scales, and Y 
is the damping rate of waves due to non-linear resonance 
with thermal and relativistic particles {Np and A^^); details 
can be found in ll24ll . 

In the following we consider a simple model of galaxy 
cluster assuming that the energy densities of relativis- 
tic protons at the beginning of reacceleration, gp, and of 
the magnetic field, B, and the injection rate of Alfven 
waves during mergers, 4, scale with thermal energy den- 
sity, £,/, (gp °^ £,/,, B and / Ii^dk °= e,/,). An exam- 
ple of the expected broad band emission (synchrotron, 
IC, n" decay) is reported in Fig. 2 by adopting the spa- 
tial distribution and physical parameters of the thermal 
IGM of the Coma cluster. In the context of this model 
the non-thermal emission is a mixture of two main spec- 
tral components: a long-living one that is emitted by sec- 



ondary particles (and by %" decay) continuously gen- 
erated during p-p collisions, and a transient component 
that is due to the re-acceleration of relativistic particles 
by MHD turbulence generated (and then dissipated) in 
cluster mergers. In order to highlight these components 
we calculate the non-thermal emission during a cluster 
merger, i.e. assuming a turbulent IGM, (Fig. 2, upper 
panels) and 1 Gyr after turbulence is dissipated (Fig. 2, 
lower panels). Lower panels show the long-living com- 
ponent of the non-thermal cluster emission, because rel- 
ativistic protons, that in turns generate secondaries, lose 
energy on a long time-scale. This long-living emission 
does not strongly depend on the dynamics of clusters 
but only on the energy content (and spectrum) of rela- 
tivistic protons in the IGM (and on the magnetic field in 
the case of the synchrotron radio emission). On the other 
hand, the comparison between upper and lower panels of 
Fig. 2 highlights the transient emission that is generated 
by short-living electrons reaccelerated by turbulence dur- 
ing cluster mergers. 

We note that the results in Fig. 2 have the potential to 
reproduce the radio bimodality observed in galaxy clus- 
ters (Fig. 1) : Radio Halos develop in connection with 
particle re-acceleration due to MHD turbulence in clus- 
ter mergers where the cluster-synchrotron emission is 
considerably boosted up (upper panel), while a fainter 
long-living radio emission from secondary electrons is 
expected to be common in clusters (lower panel); the 
level of this latter component must be consistent with 
the radio upper limits from radio observations of clusters 
with no Radio Halos. IC hard X-rays are also produced 
in connection with Radio Halos, although the IC signal 
from re-accelerated secondary electrons is not expected 
to be very luminous (see discussion in |24] for a com- 
parison with the case of re-acceleration of primary relic 
electrons). 

An important point is that gamma ray emission is ex- 
pected (at some level, depending on the content and spa- 
tial distribution of relativistic protons) to be common in 
galaxy clusters and not directly correlated with the pres- 
ence of giant Radio Halos. Cerenkov arrays already con- 
strain the level of gamma rays from Coma and few other 
nearby clusters (Fig. 2, Perkins et al., this meeting). Af- 
ter wl yr of operations FERMI will reach adequate sen- 
sitivity in the energy range 0.1-100 GeV to start obtain- 
ing crucial constraints on nearby clusters and hopefully 
to measure the energy content of relativistic protons in 
these clusters. 



ACKNOWLEDGMENTS 

We acknowledge partial support from ASI-INAF 
1/088/06/0 and PRIN-INAF2007, and useful discussions 
with M. Giroletti. 



REFERENCES 

1. Ensslin, T.A., Biermann, P.L., Klein, U., Kohle, S., 1998, 
A&A 332, 395 

2. Sarazin, C.L., 1999, ApJ 520, 529 

3. Blasi, P., 2001, Astropart. Phys. 15, 223 

4. Brunetti, G., Setti, G., Feretti, L., Giovannini, G., 2001, 
MNRAS 320, 365 

5. Brunetti, G., Blasi, P, Cassano, R., Gabici, S., 2004, 
MNRAS 350, 1174 

6. Petrosian, V., 2001, ApJ 557, 560 

7. Miniati, R, Jones, T.W., Kang, H., Ryu, D., 2001, ApJ 
562, 233 

8. Ryu, D., Kang, H., Hallman, E., Jones, T.W., 2003, ApJ 
593, 599 

9. Brunetti, G., Lazarian, A., 2007, MNRAS 378, 245 

10. Pfrommer, C, 2008, MNRAS 385, 1242 

11. Blasi, P, Gabici, S., Brunetti, G., 2007, IJMPA 22, 681 

12. Volk H.J., Aharonian F.A., Breitschwerdt D., 1996, SSRv 
75, 279 

13. Berezinsky V.S., Blasi P, Ptuskin V.S., 1997, ApJ 487, 
529 

14. Blasi, P., Colafrancesco, S., 1999, Astropart. Phys. 12, 
169 

15. Ricker, PM., Sarazin, C.L., 2001, ApJ 561, 621 

16. Dolag, K., Vazza, P., Brunetti, G., Tormen, G., 2005, 
MNRAS 364, 753 

17. Vazza, P., Tormen, G., Cassano, R., Brunetti, G., Dolag, 
K., 2006, MNRAS 369, L14 

18. lapichino, L., Niemeyer, J.C., 2008, MNRAS 388, 1089 

19. Reimer, O., Pohl, M., Sreekumar, P., Mattox, J.R., 2003, 
ApJ 588, 155 

20. Peretti, L., 2003, in Matter and Energy in Clusters of 
Galaxies, ASP Conf. Series 

21. Perrari, C, Govoni, P., Schindler, S., Bykov, A.M., 
Rephaeli, Y., 2008, SSRv 134, 93 

22. Dennison, B., 1980, ApJ 239, L93 

23. Pujita, Y, Takizawa, M., Sarazin, C.L., 2003, ApJ 584, 
190 

24. Brunetti, G., Blasi, P, 2005, MNRAS 363, 1173 

25. Giovannini, G., Tordi, M., Peretti, L., 1999, New 
Astronomy 4, 141 

26. Buote, D.A., 2001, ApJ 553, L15 

27. Cassano, R., Brunetti, G., Venturi, T, Setti, G., Dallacasa, 
D., Giacintucci, S., Bardelli, S., 2008, A&A 480, 687 

28. Venturi, T, Giacintucci, S., Dallacasa, D., Cassano, R., 
Brunetti, G., Bardelli, S., Setti, G., 2008, A&A 484, 327 

29. Brunetti, G., Venturi, T, Dallacasa, D., Cassano, R., 
Dolag, K., Giacintucci, S., Setti, G., 2007, ApJ 670, L5 

30. Dolag, K., 2006, Astron. Nachr., 327, 575 

31. Subramanian, K., Shukurov, A., Haugen, N.E.L., 2006, 
MNRAS 366, 1437 

32. Brunetti, G., Giacintucci, S., Cassano, R., Lane, W., 
Dallacasa, D., Venturi, T, Kasim, N.E., Setti, G., Cotton, 
W.D., Markevitch M., 2008, Nature 455, 944 

33. Schlickeiser, R., Sievers, A., Thiemann, H., 1987, A&A 
182, 21 

34. Thierbach, M., Klein, U., Wielebinski, R., 2003, A&A 
397, 53 

35. Cassano, R., Brunetti, G., 2005, MNRAS 357, 1313 



